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An understanding of the anomalous charge dynamics in the high-T c cuprates is obtained based 
on a model study of doped Mott insulators. The high-temperature optical conductivity is found 
to generally have a two-component structure: a Drude like part followed by a mid-infrared band. 
The scattering rate associated with the Drude part exhibits a linear-temperature dependence over 
a wide range of high temperature, while the Drude term gets progressively suppressed below a 
characteristic energy of magnetic origin as the system enters the pseudogap phase. The high-energy 
optical conductivity shows a resonancelike feature in an underdoped case and continuously evolves 
into al/u tail at higher doping, indicating that they share the same physical origin. In particular, 
such a high-energy component is closely correlated with the oj-peak structure of the density-density 
correlation function at different momenta, in systematic consistency with exact diagonalization 
results based on the t-J model. The underlying physics is attributed to the high-energy spin- 
charge separation in the model, in which the "mode coupling" responsible for the anomalous charge 
properties is not between the electrons and some collective mode but rather between new charge 
carriers, holons, and a novel topological gauge field controlled by spin dynamics, as the consequence 
of the strong short-range electron-electron Coulomb repulsion in the doped Mott insulator. 

C/2 1 PACS numbers: 74.20.Mn,72.10.-d,74.25.Gz,78.20.Bh 
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Charge dynamics in the high-T c cuprates has been under intensive studies for nearly two decades by now. With the 
O | continuous improvements in experimental techniques and sample quality, the anomalous charge dynamics exhibited 

in optical measurements has been firmly established as intrinsic properties unique to the high-T c cuprates. 
^vq | Some important features observed 1,2 in the in-plane optical measurements are: (a) The undoped parent compound 
^ . is a Mott insulator with a clear charge transfer gap observed in the optical conductivity. Upon doping, spectral 
' weight starts to appear inside the charge transfer gap. 3 (b) In heavily underdoped materials, the optical conductivity 
shows a very clear two-component feature: a low-lying Drude component followed by a resonancelike peak around the 
so-called high mid-infrared energy a>^ ir ^0.5 — 0.8 cV. 3-6 (c) Below some pseudogap temperature in the underdoped 
cuprates, the low-energy part of the optical conductivity is suppressed with decreasing temperature, accompanied by 
\^ ■ an emerging lower mid-infrared resonancelike peak at cj^ ir (~ 0.1 eV near the optimal case). 6 ' 7 (d) Near the optimal 
doping, the high mid-infrared peak around w^ ir continuously evolves into an approximate 1 /u> behavior in the optical 
conductivity, with a shift of spectral weight towards lower energy. 5,8 At the same time, the scattering rate shows a 
linear-u dependence, 9,10 while the scattering rate at w ~ increases linearly with temperature, consistent with the 
dc resistivity measurement. 11-13 

Such anomalous electromagnetic response in the cuprates has posed a great challenge to any microscopic theory 
with regard to the charge scattering mechanism. Mainly focusing on the optimal case, the one-component ap- 
proaches, including the marginal Fermi-liquid theory, 14 projected Fermi liquid with edge singularities, 15 slave-boson 
gauge theory, 16 nearly antiferromagnetic Fermi-liquid theory, 17 ' 18 etc., have all stressed how the scattering rate 1/r 
may become strongly frequency-dependent as the consequence of strong interactions. Conventional electron-boson- 
mode coupling mechanisms have been also suggested in this regard. 19-22 Alternatively the power-law behavior of the 
scattering rate \/t(oj) ~ uj a has also been interpreted as due to some zero-temperature quantum critical point hiding 
' below the superconducting dome in the phase diagram. 26 

To fit the underdoped data, however, a multi-component model composed of a free-carrier Drude term and a set 
of Lorentzian oscillators seems more appropriate based on a phenomenological consideration to account for the high 
mid-infrared band. 1 ' 2 Note that the presence of a mid-infrared band has been indeed found in various microscopic 
models describing the hole motion in a quantum antiferromagnet. 23-25 On the other hand, numerical studies 2 based 
on the t-J model and Hubbard model have also shown the existence of a mid-infrared absorption at small doping and 
its evolution into an approximately 1/uj tail with increasing doping. These all point to the fact that strong correlations 
may be essential in order to fully understand the overall charge dynamics in the cuprates. 

In this paper, we will study the global feature of charge dynamics based on an effective theory describing a doped 
Mott insulator. Due to the strong correlation nature, the charge scattering mechanism is no longer a conventional mode 
coupling between the electrons and some collective bosonic mode. It is replaced by a new type of "mode coupling" , 
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in which the spinless charge carriers, holons, interact with a topological gauge field, with the latter controlled by 
the neutral spin dynamics. Namely the strong short-range Coulomb interaction effect between the electrons is now 
represented by a novel scattering between the charge and spin degrees of freedom of the system. 

Base on this model, we show that the high-temperature optical conductivity is generically associated with a two- 
component behavior. Here the low-w Drude-like term is characterized by a linear-temperature-dependent scattering 
rate over a wide range of temperature, which gets progressively suppressed as the system enters into the pseudogap 
phase at low temperature with the emerging of a lower mid-infrared peak. On the other hand, the high-w mid-infrared 
part exhibits a resonancelike feature in the strong scattering case and continuously evolves into al/w tail in the weak 
scattering limit, which can be related to the underdoping and high doping, respectively. The high-energy mid-infrared 
component is further shown to be closely correlated with the w-peak structure of the density-density correlation 
function at different momenta which is systematically consistent with the exact diagonalization results based on the 
t-J model. Therefore, we establish a consistent picture for the charge dynamics in the cuprate superconductors in 
different temperature and energy regimes within a single unified theoretical framework. 

The rest of paper will be organized as follows. In Sec. II, the basic model describing a doped Mott insulator 
will be introduced, and various important temperature and energy scales decided by the model will be discussed. In 
particular, a novel charge scattering mechanism embedded in the model will be emphasized. In Sec. Ill, we study 
the high-temperature charge dynamics where the gauge fluctuations are reduced to static to allow for a numerical 
simulation. The optical conductivity and scattering rate are calculated. The two-component u structure and its 
correlation with the density-density correlation function will be discussed in detail. In Sec. IV, the scattering strength 
is substantially reduced in the pseudogap phase where a perturbation approach is allowed to study the low-energy 
optical conductivity. A pseudogap behavior of the optical conductivity with the emergence of a lower mid-infrared 
resonance will be determined here. Finally, a summary and discussions will be presented in Sec. V. 



From a doped-Mott-insulator point of view, the t-J model is believed to be the simplest relevant model for the 
high-Tc cuprates. Finding the correct low-energy effective theory for the t-J model has been a fascinating focus in 
search for the mechanism of high-T c superconductivity. The concepts of the resonating-valence-bond (RVB) 28 and 
spin-charge separation 29 have been conjectured in order to understand the strong correlation nature of the electrons in 
such a system, although, different from the one-dimensional case, 30 ' 31 no well-defined spinless charge carriers (holons) 
and 5=1/2 neutral spin excitations (spinons) have ever been unambiguously identified experimentally in the two- 
dimensional case. Nevertheless, there has been strong numerical evidence indicating that the charge and spin dynamics 
should be described separately by two different degrees of freedom with distinct characteristic energy scales. 30 ' 32 

Microscopic studies of the t-J model indicate that residual strong correlations still generally exist, in a spin-charge 
separation description, between the two building blocks - spinless "holons" and charge neutral "spinons" . For example, 
the interaction is mediated by the U(l) 33 or SU(2) 34 gauge field in the slave-boson approach. In contrast, in the phase 
string theory 35 the elementary force between the bosonic "holons" and "spinons" is mediated by the mutual Chern- 
Simons gauge fields. These gauge interactions complicate the issue regarding how to make a clear identification of 
holons and spinons, even if they exist in the cuprates. 

Nevertheless, the presence of a spin-charge separation will predict a set of highly nontrivial charge and spin dynamics 
whose unique features can be subjected to a systematic comparison with the experimental measurement. Even in 
low-temperature phases where the confinement of "holons" and "spinons" due to the gauge fields may occur at low 
energies, the composite structure of the electrons in terms of "holons" and "spinons" should still exhibit interesting 
features with the increase of energy, which approaches the "asymptotic freedom" at some shorter distance and higher 
energy. An experimental probe of high-energy excitations is particularly meaningful here. 

In the following, we explore the charge dynamics based on an effective gauge theory derived from the t-J model 
and show how the spin-charge separation exhibits and plays the essential role there. 



II. BASIC MODEL 



A. Phase string model 



Our starting point is the phase string model, 35 37 which takes the form: H stl i ng = Uh + H s with 




(1) 



blb)_ IJ + H.c. + \Y J [nl-{l-5)/2] 



(2) 
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in which the holon field, described by the bosonic annihilation operator hi, carries the full charge +e coupling to the 
external electromagnetic vector potential A e - . 

The spinless holon field and the neutral 5 = 1/2 spinon field (which carries the spin index as described by the 
bosonic annihilation operator 6 icr ) are minimally coupled to two internal gauge field Afj and A^-, whose gauge- invariant 
strengths are constrained to the matter fields by the following relations, respectively 

2^,. = nJ2K-n b n ) (3) 

c zes c 

c ies c 

where n\ a — bj^bia and nj 1 = h\hi, respectively, denote the on-site spinon and holon number operators, and Sc is 
the region enclosed by an arbitrary (counterclockwise) closed loop C. Such relations are known as the mutual Chern- 
Simons gauge structure which dictate that a holon and a spinon perceive each other as a fictitious 7r-flux tube [here 
to avoid the short-distance uncertainty at each center of a 7r-flux tube, on the right-hand sides (rhs) of Eqs. (3) and 
(4), the distribution of a holon or spinon at site I should be understood as being slightly smeared within a small area 
centered at /]. 

The phase string model thus defined explicitly respects all the symmetries including the translational, time-reversal, 
parity, and spin rotational symmetries. Note that the spin operators in this model are defined by Sf — 1/2 cm b il7 , 

Sf = b\^bii{— 1)V*?, and Sr = (5 , 4 + ) t , where $f — $^ = 2A , - j with the core of each flux-tube being smeared within 
a small area as mentioned above. 

In Eqs. (3) and (4), A^ and A are the chemical potentials to implement = NS and ^ i(T n\ a = N(l — $) = Nn b , 

where N is the total number of lattice sites and 6 the doping concentration. The effective hopping integral th ~ 0.67t 37 
and in this paper we shall choose th — 2 J which corresponds to t ~ 3 J. The renormalized superexchange coupling 

J s = J ( A s \ (1 — 2gS) /2 (g ~ 2) 38 with J s <~ J/2 at low doping, where t and J are the bare parameters in the original 



t-J model. Here the gauge-invariant bosonic RVB order parameter A s = \A s y = e lcrAi: >bi r7 bj- a j is determined 
sclf-consistently. 



B. Phase diagram and basic energy scales 

Fig. 1 shows the basic phase diagram 38 ' 39 for the phase string model given in Eqs. (1) and (2). The characteristic 
temperature To denotes the boundary of the so-called the upper pseudogap phase (UPP), described by the bosonic 
RVB order parameter A s ^ based on H S ; 3S The so-called lower pseudogap phase (LPP) or spontaneous vortex phase 
(SVP) is nested within the UPP with a characteristic temperature T v , which corresponds to the bosonic degenerate 
regime for the holons; 39,40 Finally the superconducting (SC) phase coherence is realized inside the LPP/SVP at lower 
temperatures. In Fig. 1 both LPP and SC phase terminate at half-filling where an antiferromagnetic (AF) long-range 
order is recovered in the ground state of H s . Note that a more careful study with considering the longer-range AF 
correlations will lead to the vanishing of the SC transition temperature T c at a finite critical doping concentration 

~ 41 

UjQ. 

Physically the UPP is a regime where short-range AF correlations start to develop quickly with decreasing 
temperature. 38 At half-filling, it continuously evolves into the AF long-range ordered state at T = 0, where the 
gapless spin-wave excitation is shown in Fig. 2 by the dotted curve. 42 However, such long-range AF correlations 
get "truncated" at finite doping after entering LPP, where an energy gap E g is opened up around the AF wave 
vector Q = (it, it) as shown in Fig. 2, which depicts the peak positions of the spin dynamic susceptibility function 
Imx zz (q, lo) at T = 42 (see Appendix A). 

The doping dependence of the resonancelike spin energy E g is illustrated in Fig. 1 which scales with T c as T c ~ 
Eg/Ak-Q as previously established 43 based on the phase string model. The peak structure of Imx zz (q, u>) in Fig. 2 
coincides with the spin- wave spectrum at half-filing (dotted curve) and becomes a non-propagating mode with a finite, 
doping-dependent width in the momentum space at finite doping. But the overall high-energy envelop of the latter 
still well tracks that of the (softened) spin wave, implying that the short-range AF correlations remain strong in the 
short-distance, high-energy regime in the UPP. This is consistent with the neutron scattering measurements. 44 There 
is a characteristic upper-bound energy scale, E^ pper , indicated by the arrows in Fig. 2 for the S = 1 excitations. 
Its doping dependence follows closely with To in Fig. 1 and monotonically increases with reducing doping, reaching 
£u PP cr _ 2.3J at half-filling. 

The low- lying sharp resonancelike structure in Inr\; zz (q, ui) at E g is the consequence 42 of the holon condensation 
realized in Hh which then influences Eq. (2) via the gauge field A\, in terms of Eq. (4). With the increase of 
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FIG. 1: In the phase diagram of the phase string model, the characteristic temperature To for the upper pseudogap phase 
(UPP) and T c for the superconducting phase (SC) are shown to be well scaled with two basic energy scales, _E" pper and E g (see 
Fig. 2 for their definitions), with the proportional coefficients ~ 0.4 and 0.25, respectively. The lower pseudogap phase (LPP) 
or the spontaneous vortex phase (SVP) is characterized by a temperature T v which is also determined by the spin excitation 

39 

spectrum. 

temperature, the holon condensation (or more precisely the amplitude condensation of the holon field) will terminate, 
not at T c , but at an intermediate characteristic temperature scale T v lying between T and T c that represents the 
onset of the LPP/SVP in Fig. 1. Such a bosonic degenerate regime corresponds to the Nernst regime observed in the 
cuprates. 45 



The charge dynamics is governed by Hh [Eq. (1)], in which a bosonic holon, carrying the full charge +e, will 
get scattered off by an internal gauge field A\y According to Eq. (3), an isolated spinon excitation will serve as a 
7r-flux tube perceived by the holons and thus provides a strong, unconventional charge scattering source. So at high 
temperature when a lot of spinons are thermally excited, one expects a severe intrinsic frustration effect exerted from 
A*j on the holons. On the other hand, the effect of the 7r-flux tubes bound to those spinons which are RVB paired at 
short-distance will be essentially cancelled out. At sufficiently low temperature and low energy where the majority of 
spins remain short-range RVB paired, the gauge fluctuations in A\- will then be substantially reduced to result in a 
weak scattering, which warrants a perturbation treatment. 

To sec how the spin dynamics influences the charge degree of freedom via Afj, one may introduce A s • (r — r^) = A?- 
and express the propagator of the gauge field A s in the continuum limit as follows 



using Eq. (3), where x zz (q, ^n) is the z-component spin susceptibility function. 

The detailed spin dynamics described by H s in Eq. (2) has been systematically studied 38,42 before. For example, 
the peak structure in Imx Z2 (q, uu) at the mean-field level 42 (T = 0) is presented in Fig. 2 at various dopings. At finite 



C. Novel scattering mechanism 




(5) 
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FIG. 2: Basic energy scales, E™ ppBT and E g , determined by the dynamic spin susceptibility function Imx zz (q, w) at T = 0. The 
peak positions of Imx zz at 8 — is shown in the energy and momentum (along the q x — q y axis) by the dotted curve, which 
tracks the spin wave dispersion with £'" ppcr ~ 2.3 J and E g = 0. The upper-bound energy _E" ppcr monotonically decreases with 
increasing doping from 0.05, 0.125, to 0.2. E g denotes the resonancelike peak energy at q = (n, n), which emerges in the LPP 
and SC state at finite doping. Note that the finite horizontal bars at finite doping indicate the momentum widths for these 
non-propagating modes. 42 



doping, due to the opening up of the low-lying spin gap E g , one finds a vanishing spectral weight of ImD A (q, uj) 
at uj < Eg and T — according to Eq. (5). Since the spin gap E g is directly related to the longest-size RVB pairs 
of spins, the effect of A s , contributed by the it flux tubes bound to individual spinons in terms of Eq. (3), will get 
cancelled out cleanly in such a low-energy, long- wavelength regime. Consequently the bosonic holons in Eq. (1) will 
be free from the flux frustration and thus experience a Bose condensation at a sufficiently low temperature, which 
defines T v for the LPP as previously stated. In fact, the zero-temperature spin dynamic susceptibility in Fig. 2 is 
obtained based on H a under the holon condensation. Namely, the spin gap and holon condensation actually enforce 
each other self-consistently in the phase string model. 

The gauge fluctuations in fl A ' (q, uj) will gain a finite strength at E g < uj < E^ ppcr and then fall off rapidly 
beyond E^ ppcr based on Imx Z2 (q, oj). At w ^ Ef ppcr , the holons will mainly experience quasistatic gauge fluctuations 
concentrated at u> < £'" ppcr , with the total strength determined by 

(<S>hf)= [ ^1^47r 2 ^(q, W ) (6) 
J q 

where $□ = a 2 z- (V x A s ) defines the local flux per plaquette (surrounding a lattice site) and the spin structure 
factor 5 22 (q, uj) = n^ 1 [1 + n (uj)] Imx zz (q, u>). According to the sum rule discussed in Ref. 42 , one finds 



(*>)*)* JQ-W-Q^V] (7) 

V □/ / y 3 I pi max v ' 

In particular, if the temperature is further increased to T > T , there is no more significant AF correlations among 
spins as A s = such that (see Appendix A) 

S zz ( q ,Lj) = -± 1 \<S>°\ 2 6(u>) (8) 

v ' 47T max ^ ' v ' 



Then the corresponding gauge flux fluctuation becomes truly static with the weight y\(^p) / concen trating at 
oj = 0. 

Namely, the gauge field A s simply describes the randomly distributed static flux at T > To and quasistatic flux at 
high energy ui ^> E" ppcr in the case of T < To. Such regimes represent the maximal quantum frustration that the 
holons can experience in the phase string model governed by Eq. (1). 
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III. CHARGE DYNAMICS AT HIGH TEMPERATURE: TWO-COMPONENT FEATURE 

In Sec. II C, we have seen that the strongest scattering to the charge carriers in the phase string model will set in 
at T > T above the UPP, where the fluctuations of the gauge field A s are concentrated at uj — 0, i.e., in the static 
limit. In the following we shall study the charge response as characterized by the optical conductivity and density- 
density correlation function in this regime. It has been noted that even at T < T , if w > £7 upper , a quasi-static 
approximation for the gauge field Afj should also be justified. In other words, the results obtained at T > T should 
be still qualitatively applicable at T < Tq at sufficiently high energies. 



A. The density of states for holons 

In the phase string model, the charge holons are coupled to the spin degrees of freedom via the topological gauge 
field Afj defined in Eq. (3). Thus the charge dynamics is strongly influenced by the spin dynamics which decides the 
propagator of Afj in terms of Eq. (5). As discussed in Sec. II C, the spin fluctuations are greatly reduced in the 
UPP at T > T and become "static" [Eq. (8)] at the mean-field level as A s = 0. Accordingly, A\- describes a static 
random flux at T > T with the strength given by Eqs. (6) and (8). 

Namely, above the UPP, the holons in Hh simply see a collection of ±7r flux-tubes bound to spinons of number 
(1—6)N which are de-paired (A s = 0) and randomly distributed in this regime with neglecting the spin-spin correlation 
at the mean- field level. 

Then, at T > Tq one can make a direct numerical simulation to diagonalize Hh for each random configuration of 
spinons like an impurity problem by a unitary transformation hi — ^ m Ci m a m with Ci m satisfying 

-t h ]T e iA °iC jm = e h m C im (9) 

j—N N (i) 

where j — NN(i) denotes the four nearest neighbor sites. Note that the hard-core interaction between the holons is 
neglected here as we are mainly interested in the high-temperature or high-energy charge behavior below. 

For the simplicity of numerical simulation, we define the smeared flux strength of Af • on each lattice plaquette by 

E A h=*n (10) 

plaquette 



and assume $ n = ± |$ n | max with the signs + and — randomly distributed on the lattice plaquette and 

|*J = (1 - <5)tt (11) 

I L-l I max v ' v y 

Note that this assumption is valid when the individual 7r-flux tubes of total number of 1 — S spinons are well dis- 
tinguished at high energy, short-distance scales, and |$ n | max is compatible in magnitude with |$Q| m in Eq- (7) 
of the mean-field version. But it should be emphasized that |$ n | max in Eq. (11) only provides an upper bound for 
the average flux strength per plaquette. One may think of many effects that can reduce its magnitude. For example, 
two flux tubes of opposite signs sitting at the nearest sites can partially cancel each other; the residual weak AF spin 
correlations among spinons will further enhance such cancellation; in addition, at larger doping the hopping of holons 
is always accompanied by a "backflow" of spinons which can mix ±n fluxes and reduce the average strength of |$ n | 
from I I 

I u I max 

Fig. 3 illustrates the holon density of states (DOS) obtained at S = 0.125 (solid curve). Here we have used the 
quenched method to average static random flux configurations of $ n . The result shows how the DOS drastically 
reshaped by the gauge field, i.e., the suppression in the high-energy (mid-band) DOS, as compared to the flux free 
case (dashed curve). Note that the dotted curve in Fig. 3 represents the DOS for the case of a uniform n flux per 
plaquette, which looks similar to the present random flux case except that the momenta remains well defined in a 
reduced Brillouin zone in contrast to a strong mixing of momenta over a wide range by the scattering effect in the 
latter. By comparison, the DOS with a reduced |^ n | m = OAtt (dash-dotted curve) is also presented. It is noted 
that other kinds of treatment for the static gauge field Afj will generically lead to the similar overall behavior. For 
example, one may treat 3> n as a white- noise random flux distributed within the interval [— |$ n | max , + |<& n | max ] or 
even introduce some spatial correlation based on Eq. (5). 
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FIG. 3: Solid curve: the holon density of states (DOS) determined by //^with the flux strength decided by Eq. (11) at 
8 = 0.125 (see text) and the dash-dotted curve: with the reduced flux strength |$ n | = 0.4-7T. By comparison, the dashed curve 
represents the flux- free limit, while the dotted curve corresponds to the case in the presence of uniform 7r flux per plaquette. 
The calculation is done on a 32 x 32 lattice. 



B. Optical conductivity: Two-component structure 

Let us consider the optical conductivity using the Kubo formula 

<7 xx (q = 0, w) = - [n M (q = 0, w) - e 2 (T x )] (12) 

where IF^q = 0,u) = U xx (iuj n — > u + i0 + ) is the retarded current-current correlation function with H xx (iuj n ) = 
— dTe tuJnT (T t J x {t) J x (0)} and the charge current operator defined by 

J x (q) = ^§£ e " q ' r * - ftj^e^.^) (13) 



(T x ) in Eq. (12) is defined by 



(14) 



The real part of the optical conductivity is then given by 

<4» = £ £ m toto , ~ n( ^' } 5(w - z m > + U (is) 



m,m' >m ^ 

in which n(£ m ) = l/(e /3 ^ m — 1) is the Bose distribution factor with £ TO = e' r l n + Xh and the matrix element 



M mm , = eH\ 



'^{e bA ^- i C* +x m C itm i - e %Ai - i +" C* m C i+XtTn >) 



(16) 



with using Eq. (9) for a given static configuration of A\-. Note that the final <j' xx (lj) is obtained by averaging over 
the quenched random flux configurations as discussed above. 

Fig. 4 shows the calculated real part of the optical conductivity at 6 = 0.065, with choosing the maximal flux 
strength in terms of Eq. (11). The main feature of the spectral curves at various temperatures T > T is that they 
can all be decomposed into a two-component structure with a usual low-energy Drude component (~ l/^ 2 ) and a 
mid-infrared resonancelike peak around the energy scale w^ ir ~ 4th, with the whole spectrum eventually terminated 
below the energy Qth- 
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FIG. 4: The real part of the optical conductivity at 5 — 0.0625 with ]$□ | = |<3?n | max and To ~ 0.5th- 




FIG. 5: The real part of the optical conductivity at 8 — 0.125 with \&n \ chosen from |$p| max = 0.875-7T to 0.2-7T at a fixed 
T = 0M h - 



1. Mid-infrared peak 

The origin of the mid-infrared resonance has been one of the most intriguing optical properties in the underdoped 
cuprates. 1 ' 2 ' 5 Normally a photon with the momentum q ~ cannot excite a high-energy particle-hole pair involving a 
large momentum transfer due to the momentum conservation law, unless there is a scattering mechanism to strongly 
and widely smear the momentum. This is difficult to realize in a conventional clectron-collectivc-mode coupling 
mechanism. The present model provides an alternative scattering mechanism due to the strong correlation effect 
caused by the on-site Coulomb repulsion in a doped Mott insulator. 

We have already seen that the effect of A s results in a double-peak structure in the holon DOS (Fig. 3). In contrast 
to the uniform 7r flux case shown in the same figure, which also has a double-peak structure, the high-energy inter-peak 
transition at q — > becomes possible in the random flux case due to the mixing between the small and large momenta 
by the strong scattering via A s . This is the origin for the mid-infrared peak found in Fig. 4. 

The presence of a peak in the optical conductivity around to ~ 3t — At has been previously identified in the Hubbard 
and t-J models by exact numerical simulations 27 ' 46 Two approaches are consistent in this regard. Note that for such 
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FIG. 6: The scattering rate l/r(w) defined by Eq. (17) at various temperatures between To ~ 0.25th and 2To which show a 
rough linear-o; dependence over a wide range at to > feeTo- Inset: 1/t(lo) vs. ui at different |$rj|'s corresponding to Fig. 5. 

a high-energy and short-distance physics the finite-size effect in the exact diagonalization calculation should not be 
important. In Sec. Ill C below, we shall further compare the density-density correlation function obtained by the 
exact diagonalization and present approach, where the w-structures at different momenta can provide much richer 
features in further support of the consistency. 

The above-discussed mid-infrared peak in the optical conductivity seems in contrast with the approximate l/u> 
behavior observed experimentally for the optimally and over-doped cuprates. 5 But if one artificially reduces the 
magnitude |$ n | from given in Eq. (11), the mid-infrared peak will actually smoothly evolve into the 1/ui 

behavior with reducing $ p |, as clearly illustrated in Fig. 5 at a fixed holon concentration S = 0.125. In Fig. 5, the 
mid-infrared resonancelike peak at smaller | $ n | 's becomes softened and finally behaves like al/u tail in the regime 
~ 2th — &th with the weight shifting towards the lower energy. As explained before, the strength of |$ n | is expected 
to decrease faster with increasing doping than |$ n | max given in Eq. (11). Such a picture is also consistent with the 
exact diagonalization calculations in the t-J model. 47 



Now let us closely examine the low-w component of the optical conductivity. Experimentally the scattering rate is 
normally defined by 



which is determined by the measured optical conductivity. Here lo p denotes the plasma frequency, which in the present 
case is given by lo p = y / 8i:e 2 dtf l . 

In Fig. 6, 1/t(w) based on the calculated a(uj) is plotted as a function of uo in different temperatures at S = 0.125. 
Here we have chosen |$ n | = 0.4-7T which corresponds to the case where the high-a; optical conductivity looks more 
like a 1/u behavior (Fig. 5). Here one finds that 1/t(u) increases monotonically with u> and is roughly linear-w 
dependent over a wide u> region at to > k^To. Note that generally the w-dependence of 1/t(u>) at higher energies 
is closely correlated with the evolution of the aforementioned mid-infrared feature, as shown in the inset of Fig. 6. 
Especially, as shown in the inset of Fig. 6, the non-monotonic behavior of l/r(w) is also found in the strong gauge 
fluctuation case (e.g., |$ n | = 0.8757r), which corresponds to the heavily underdoped case, and is consistent with 
recent experimental results. 2,5 

In particular, one sees a parallel shift of 1/t(u>) with increasing temperature at 1ow-cj in Fig. 6, which implies a 
linear-temperature dependence of the dc scattering rate. Such a parallel shift is also observed in the experiments 10 ' 
which persists to a very high energy scale (~ 3000cm _1 ), implying an unconventional scattering mechanism in such 



2. Low-energy component 




(17) 
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FIG. 7: The dc resistivity p dc — l/cr'(w)|o>^o as a function of temperature which is fit by a straight line showing the linear-T 
dependence. Inset: p dc at different |<&rj|'s which all show good linear-T behavior with slightly different slopes. 



strongly correlated systems. The dc scattering rate 1/Vd c can be determined by extrapolating 1/t(oj) to ui = 0. Due 
to the parallel shift, the temperature dependence of 1/t(o>) at low-u is quite similar. The obtained dc resistivity 
based on the Drude formula p dc = (uj^/Ah) r^ 1 = l/er'(0) is shown in Fig. 7 which is indeed quite linear over a very 
wide range of temperature at T > T . 

In order to understand the physical origin of the linear-T behavior, one may rewrite cr'(0) in terms of Eq. (15) as 
follows 



^'(0)=/^£MO[l + n(CJ] 



where the chemical potential in £ m is determined by n(£ m ) = NS. We find cr'(0) oc (5 over a very wide range of 
the temperature at T > T where n(£ m ) -C 1, i.e., in the classical regime of the bosons. The corresponding scattering 
rate h/r& c ~ 0.7fceT for the case shown in the main panel of Fig. 7, whose slope is slightly |<fr n | dependent as 
indicated in the inset. Indeed, as discussed in Ref. 39 , the bosonic degenerate regime already ends up at T v , i.e., the 
boundary of the LPP/SVP (Fig. 1), where the excited spinon number becomes equal to the holon number such that 
the quantum phase coherence among the latter get totally interrupted by the former which carry 7r-fiux tubes. At 
T > To, totally 1 — 8 randomly distributed 7r-nux tubes are perceived by S holons and the latter behave like classical 
particles. One expects this anomalous transport be smoothly connected to the Brinkman-Rice retracing path regime 49 
in the large T limit. 

The dc scattering rate fi/rd c ~ 2fc^T has been previously obtained 50 by the quantum Monte Carlo numerical 
method, where the starting model is a system of interacting bosons coupled with strong Gaussian fluctuations of the 

static gauge field of the strength ^(^p) 2 ^) • Note that ^(^q) 2 ^) used in the Monte Carlo simulation is about the same 
order of magnitude as in the above case and, in particular, it is temperature independent in contrast to the linear-T 
dependence predicted in the slave-boson U(l) gauge theory 16 which was the original motivation for such a Monte 
Carlo study. 50 It is noted that the short-range repulsion between the holons has been taken into account in the Mont 
Carlo calculation, which may be responsible for the larger temperature slope of the scattering rate as compared to 
the present approach besides the Gaussian approximation used there. A further difference in the treatment is that 
we have used a quenched method to average over the static random flux configurations of $ n since Afj depicts 7r-flux 
tubes bound to random distributed spinons at T > To, while an annealing approximation is used in Ref. 50 because of 
the Gaussian fluctuations of the gauge flux employed there. As we are mainly interested in the high-temperature (or 
high-energy) behavior, such a difference is qualitatively not important as discussed in Ref. 51 . 

It is noted that the lattice effect becomes very important at high-energy, short-distance at such strong flux fluctu- 
ations. It is actually responsible for the double-peak holon DOS, the mid-infrared feature in the optical conductivity, 
and the high-w behavior of 1/t(u>) discussed above. In the following we examine an another consequence of this 
unique DOS structure. 
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FIG. 8: (a) The density-density correlation function Cd(fl,(jj)/8n at q = (tt,tt/3) at S — 0.0625 at various temperatures; (b) 
The scaling behavior of C<j(q, u})/8n at different dopings at T = 0.5th- 



C. Density-density correlation function 



In the above section, the mid-infrared resonance peak of the q = optical conductivity has been attributed to a 
large-w transition between the double peaks of the holon DOS [Fig. 3], which is the consequence of the holon coupling 
with the strong fluctuating gauge field. In the following we discuss an independent probe of such a peculiar DOS 
structure by studying the density-density function at finite momentum q and energy to, and compare the results with 
the exact numerical calculations. 

The charge (holon) number operator in the momentum space is expressed as 



i{c X ) = Y,^ ri h\h i 



(18) 



Similar to the (retarded) current-current correlation function, the imaginary part of the (retarded) density-density 
correlation function can be expressed as 



m,m' 



K£m) - «(£m')] ~ Cm' + O 



Corresponding to the calculated optical conductivity in Figs. 4 and 5, using the same parameters and method, the 
structure function Cd(q, a;) of the density-density correlation can be similarly computed. 

The calculated Cd(q, lu)/5tt at momentum q = (w, 7r/3) are presented in Fig. 8 at different temperatures and doping 
concentrations. The double-peak structure as a function of u> reflects the underlying structure of the DOS of holons 
in Fig. 3 with $ n = ± <& p An interesting feature shown in Fig. 8(b) is that C^q, iv)/5ir is roughly doping- 

independent. Such a "scaling" behavior was previously found numerically for the t-J model, 32 which is generally 
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FIG. 9: The calculated density-density correlation function Cd(q, u))/5t[ at different momenta (solid curves) with T = Q.bth- 
The exact diagonalization results 32 are shown as dashed curves. 



inconsistent with the picture of particle-hole excitations in a Fermi liquid and was conjectured to be an indication of 
the bosonic description of charge excitations. 32 

The calculated Cd(q, u))/5ir is presented in Fig. 9 (solid curves), which evolves distinctively with different momenta. 
For comparison, the exact diagonalization results 32 are presented as dotted curves. Due to the approximate doping- 
indcpcndence, the numerical result of four holes in a 4 x 4 lattice is used here. It is interesting to see that the overall 
w-peak feature of the calculated density-density correlation function is in qualitative and systematic agreement with 
the numerical one at different q's without fitting parameters (here t is simply set at th as the mid-infrared feature 
peaks around ~ At in the numerical calculation). Such a consistency between the present effective theory and the 
exact diagonalization provides an another strong evidence, in addition to the mid-infrared resonance peak of the same 
origin, that the gauge-coupling boson model (1) correctly captures the high-energy charge excitations in the t-J model 
and large- U Hubbard model. 



IV. LOW-ENERGY PSEUDOGAP BEHAVIOR AT LOW TEMPERATURE 



So far we have been focused on the "normal state" above the UPP at T > Tq where we have seen the maximal 
static scattering coming from the gauge field. As outlined in Sec. II, when the temperature decreases below To, the 
spins start to form bosonic RVB pairs and develop short-range AF correlations. Although this does not change the 
overall integrated strength of the gauge fluctuations, due to the spin sum rule as given in Eqs. (6) and (7), the low-w 
gauge fluctuations of A s will get suppressed progressively. Namely the scattering felt by the low-lying charge carriers 
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q,ip„ 



kjco k'.ito ' k.ito 



FIG. 10: Feynman diagram for the boson self-energy S(k, iio n ). 

will be reduced. Here the gauge fluctuations are no longer static and are strongly correlated with the spin dynamics, 
where the above static approximation method is not applicable at least in the low-energy regime. 

In the following we shall study the low-energy optical conductivity at T < T by using the perturbative method 
and continuum approximation instead. Such an approach is meaningful in the regime where the spin fluctuations 
are substantially suppressed, which in turn results in the weak fluctuations of A s according to Eq. (5). Of course, 
as emphasized before, in the high-energy regime (u> 3> Ef pper ) the charge carriers still feel the strong scattering by 
a quasistatic gauge fluctuation of the overall strength given by Eq. (7) and the previously results at T > T are 
expected to still hold qualitatively. 

A. Self-energy 

In the case of weak gauge fluctuations, one can take the continuum limit in the Hamiltonian (1) as follows 

H h = -^J d 2 vh\v) [V - iK s (v)f h(r) + X h j d 2 vh\v)h(v) (19) 
where rrih = (2i/ja 2 ) 1 with a shift in the chemical potential A/j — Atu — ► A^. Further express Hh = H a + H i with 



* - 1 



2m h 



t&h v (21) 



in which q = k k' and qi + q2 = k — k'. 

We treat Hi as a perturbation and calculate the holon self-energy up to the quadratic order of the gauge field. In 
the imaginary-time representation, the self-energy can be written as 

£ (k ,^„) = --^_|] / ^G( )(k',zc; i )^(q, i p„)(k + k'r(k + k')' 3 (22) 

where the free boson field propagator G^(k, iuj n ) = 1/ (iuj n — £ k ) with £ k = + X h . The corresponding Feynman 
diagram is illustrated in Fig. 10 with q = k k' and p n = u n — ui' n . By using the spectral function representation 
of D A in terms of lmx zz according to Eq. (5) and after the frequency summation and analytic continuation iuo n — > 
lu + i0 + , one finally obtains the imaginary part of the self-energy 

Att 2 f°° f r/ 2 k' Ik x nl 2 

S"(k, u) = ~J^2 J *»' 'J J^ lm * ZZ ^ w ')^i 5L [»KkO + n(u>') + 1] 5(cj £ v J) (23) 

Based on the dynamic spin susceptibility function lmx zz determined 42 by H s (see Appendix A), E"(k,£ k ) can be 
then numerically computed and the results are presented in Fig. 11 at several low-temperatures below T . It shows 
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FIG. 11: The imaginary holon self-energy at various temperatures below Tq. 



that generally £"(k, £ k ) oc £ k at high energy with the slope dependent on temperature. But at low energy, E"(k, £ k ) 
is quickly suppressed below Ef ppcr at T = or reduced from the linear-e k behavior at finite temperature due to the 
presence of spin dynamics, characterized by lmx zz shown in Fig. 2, which makes the total strength of the gauge 
fluctuations spread over a finite w region such that the effective scattering becomes weakened at low energy. 

However, in contrast to a fermion system, the self-energy E"(k, £ k ) is not generally related to the measurable 
transport properties for the present boson system. So in the follow we shall directly consider the optical conductivity 
based on the Kubo formula. 

B. Optical conductivity 

The optical conductivity <j'(u>) is determined by the retarded current-current correlation function. Note that in the 
continuum limit, the gauge-invariant current in the phase string model is given by 



which can be further written in two parts J = J A + J B , with 

e 



J d 2 r[tf{v)Vh{v) - iA s {r)h(r)^ h{r)] (24) 



ja = ^y kh t hk (25) 

JB = -^P'^ (26) 

K.q 

Correspondingly, to leading order of approximation, the current-current correlation function contains two parts: 
n = n A + n B , where 

e 2 

U A (iu n ) = - m 2 Na 2 fc2 12 Gh ^ ipn + iuj n)G h (k, i Pn ) (27) 
nB (^«) = E D A \c l MG%\k,ip n -iu; n )G£\k + c l ,ip n + ip' n ) (28) 

n ' k,q iPn,iP n 

with G(k, iu n ) = l/ [iu n - £ k - E(k, iiv n )] . 
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FIG. 12: Optical conductivity at different temperatures below To. A new lower mid-infrared peak emerging at w^ir which is 
weighted between the magnetic energy scales E g and _E" ppcr . 



After the frequency summation and analytic continuation, one finds 

M A (.) = -^TJ^ E fc2 J T* Ah[K £)Ah[K " + a ' )[n(£) _ n{£ + W)] 



h k 

where the spectral function Ah (k, is defined as 

-2E"(k,w) 



A(k,w) 



[c-£ k -E'(k,^)] 2 + [E"(k, W )] 2 



It is easy to find that without including E(k, w), ImIl A (u)) = at finite lo. The behavior of E"(k, uo) has been 
discussed in the last subsection. Its correction to Imn j4 (cj) turns out to be in a higher-order as compared to the 
leading contribution in ImW B {lo) given below 

ImIT» = Yl / du'hnD A '(q,u')[n(^) n(£ k+q )] {[n( W - u/) + n(a/) + 1] 

x % - (£ k+q - £ k ) - u/) + [n(w') - n(w + w')]<f (w - (£ k+q - 4) + w') } (29) 



To study the low-w behavior, thus one needs only to keep the second term: cr'(u)) w — ImIl B (ui) / uj . In Fig. 12, 
we have plotted the calculated optical conductivity at various low temperature below To at S = 0.125. A prominent 
suppression of a'(ui) at low-w is present at T = with a second "mid-infrared" peak emerging around w^ ir ~ 0.75 J 
which sits somewhat between the two characteristic magnetic energy scales, E g and _E" ppcr , as marked in the figure. 
Note that such a new energy scale in the low-w optical conductivity merely reflects some weighted energy scale based 
on the magnetic Imx zz , which may be seen from the following simplified formula for the optical conductivity at low 
temperature 

With the increase of temperature, the "gap" at low energy in cr'(w) is quickly filled up by the thermal excitations as 
shown in Fig. 12. The lower "mid-infrared" peak feature remains around u>^ ir at low temperature throughout the 
LPP below T v . Note that T v is between T c and To, and the dashed curve at T = 0.4T is obtained by supposing that 
T > T v where Im\ zz behaves differently. 38 As compared to the solid curve at the same T = 0.4T , which corresponds 
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FIG. 13: A schematic optical conductivity summarizing the general behavior at different energy and temperature regimes in 
the phase string model. 



to the case inside the LPP, the overall difference is small except for the vanishing the lower "mid-infrared" peak [Fig. 



Finally, we note that a Drude behavior without any signature of pseudogap feature is indeed recovered if the phase 
at T > To is extrapolated to T = with using Eq. (8): 



Of course, the perturbative approach is no longer expected to work reliably at such strong flux fluctuations. 



We have studied the charge dynamics based on the phase string model of the doped Mott insulator. This model 
possesses a mutual Chern-Simons gauge interaction between the bosonic holons and spinons, which is responsible for 
some highly non-trivial charge and spin dynamics in the system. 

In the present work, the optical conductivity generally exhibits a two-component feature: a "coherent" low-energy 
part plus an "incoherent" high-energy part. The former is Drude-like at high temperature, which gradually gets 
suppressed, after entering the upper and lower pseudogap phases, at to < w^ ir . This is related to the reduction of 
spin fluctuations, with a lower "mid-infrared" peak emerging at in o~'(u>) which is originated from a weighted 
magnetic energy scale between E g (associated with the resonancelike energy in the dynamic spin susceptibility) and 
^uppcr ^ an U pp er bound magnetic energy for S = 1 excitations). The high-energy incoherent part shows a resonance 
around w^ ir ~ 4th if the gauge fluctuation is sufficiently strong, which can smoothly evolve into al/w tail at weak 
gauge fluctuations, consistent with the situations in undcrdoping and optimal/over-doping, respectively. The overall 
features of the optical conductivity are summarized in Fig. 13, whose behavior is specified by two energy scales: w^ir 
and oj^ ir , which have been observed 7 experimentally. We have also shown that the high mid- infrared resonance at w^ ir 
can be closely correlated with the ui structure in the density-density correlation function which agrees qualitatively 
with the exact diagonalization results at all momenta. 

The high-temperature Drude component is characterized by a linear-T scattering rate l/r^c in the region of T > Tq. 
This is consistent with the previous result obtained 50 using the quantum Monte Carlo simulation for a phenomeno- 
logical model of repulsive bosons interacting with a strong spatially fluctuating gauge field. Furthermore, in the same 
high-temperature regime, the scattering rate 1/t(u>) also shows an approximate linear to dependence in the case when 
the high-w optical conductivity evolves into al/u behavior. Thus, the phase string model can give a consistent expla- 
nation for the anomalous dc and ac transport in the normal state of the optimally doped high-T c cuprates. Note that 



12]. 




(30) 
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such linear Toru dependence of the scattering rate above the pseudogap phase is really a high temperature/energy 
behavior and thus is distinct from a quantum critical phenomenon. 

Physically, the phase boundary at To separates the high-temperature classical regime from the spin "degenerate 
regime" in the UPP where the spins form the RVB pairing with A s ^ 0. In the classical regime, both the charge 
and spin degrees of freedom behave diffusively, and the linear-T relaxation rate may be regarded as the system in the 
"quantum limit of dissipation" . 52 It is interesting to see how small To is as compared to a normal Fermi degenerate 
temperature scale. The charge holons will gain quantum coherence below a lower temperature T v in the "boson 
degenerate regime" as illustrated in Fig. 1, where the optical conductivity exhibits the pseudogap behavior at low 
energy. Note that at T > To, although A s = 0, the superexchange J still plays a critical role 35 to "repair" the 
spin "z-componcnt mismatch" created by the hopping of the holes to ensure an effective hopping integral, while the 
"irreparable" transverse spin mismatch described by the phase string effect is responsible for the maximal gauge flux 
frustration discussed in this paper. Such a "normal state" is expected to be smoothly connected to the Brinkman-Ricc 
retracable path limit at T > To, where the role of J eventually diminishes and the present effective model is no 
longer valid. 

Therefore, the simple boson model (1) seems to well capture many important high-energy properties of the charge 
excitations in the cuprates and the t — J model. Another interesting channel to probe the high-energy spin-charge 
separation is the ARPES experiment. We expect to see the distinct energy scales in the single electron spectral function 
and a detailed study of the high-energy structure based on the phase string model will be presented elsewhere. It is 
further noted that in the present work we have not considered the quasiparticlc contribution to the optical conductivity. 
In the phase string model, a spin-charge recombination will occur in the superconducting phase in which nodal 
quasiparticles become stable below the energy scale set by E g . 53 But this is a fairly low energy as compared to to^ ir , 
which should not change the overall picture presented in Fig. 13. 



APPENDIX A: SPIN DYNAMIC SUSCEPTIBILITY FUNCTION 

According to Ref. 42 , the spinon Hamiltonian H s in Eq. (2) can be diagonalized as H s = J2ma E m lmalma + const, 
by the Bogoliubov transformation 

ha = w rna(i)[u,n1 m a ~ "mll-J (Al) 

m 

with the spinon wave function uj ma (i) determined by the eigen equation: 

£ m uw(i) = -J s £ e-^uw(j) (A2) 

j=nn(i) 

and 



""' = 7lVit + 1 <A3) 



^HWib- 1 (A4) 



where the spinon excitation spectrum E m = y A 2 — £ m . Here the Lagrangian multiplier A is determined by enforcing 
the average constraint (j^a ^L^) = 1 — ^> leading to 

m 

Correspondingly the dynamic spin susceptibility function can be expressed as 42 

r0 



X-(q,iu,„) = ^E/ {StWSfie-^-'A+^dT 
~7 ^ ' C mm i a (q) 



4 

mm' a 



, _ ^2 n{E m ,) - n(E m ) 

lUJ n + tL m — ±b m i 



+ (u m v m > - v m u ml ) (1 + n(E m ) + n{E m ,)) x - 



2 V iui n + E m + E m ' iu) n — E m — E m i 



(A6) 
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where C mm > a (q) = jfj^ij e * q ' (r * rj)uj * ma (i) w m« {i)^m'a (O w m'a CO- By taking the analytic continuation iu n -> 
u> + z0 + , one finally obtains 

Imx Z2: (q,w) = ^ C mm , a (q) {(u m u m / - v m v m ,f [n{E m >) - n{E m )] 8 {w + E m - E m >) 

mm' a. 

2 1 

+ (u m v m > - v m u ml ) [1 + n(E m ) + n(E m >)\ x - [5 (w - £ TO - £ m -) - (5 (w + £ m + E rn ,)]} 

In the UPP at T > T , J s — > such that _E m = A, w m = 1 and w m = 0. A is determined by Eq. (A5) as 38 
e /3A = (3 — J) / (i — J). i n this case, one finds the spin structure factor 

S"(d,w) = 7r- 1 [l + n(o;)]Inix"(q,w) 

= ^(3-i)/(l-i)%) (A7) 

which is concentrated at w = as there is no spin-spin correlation at the mean-field level. Note that in obtaining the 
last line a correction factor 2/3 is multiplied such that the precise sum rule J dw J^q S zz (q_, oj) — N (S z ) 2 — N/4 is 
satisfied at half-filling. 42 

On the other hand, at T = 0, due to the holon condensation, A'^ can be treated as describing a uniform flux 
satisfying 

E A % = ^ ( A8 ) 

plaquette 

which is expected to persist up to T v . The detailed solution has been presented in Rcf. 42 . In Fig. 2, the peak positions 
in Imx zz (q, oj) at different doping concentrations are shown. 
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